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The phenomenological classification of collective quadrupole excitations by means of 
the Bohr Hamiltonian is reviewed with focus on signatures for triaxility. The variants 
of the microscopic Bohr Hamiltonian derived by means of the Adiabatic Time Depen¬ 
dent Mean Field theory from the Pairing plus Quadrupole-Quadrupole interaction, the 
Shell Correction Method, the Skyrme Energy Density Functional, the Relativistic Mean 
Field Theory, and the Gogny interaction are discussed and applications to concrete nu¬ 
clides reviewed. The Generator Coordinate Method for the five dimensional quadrupole 
deformation space and first applications to triaxial nuclei are presented. The phenomeno¬ 
logical classification in the framework of the Interacting Boson Model is discussed with 
a critical view on the boson number counting rule. The recent success in calculating the 
model parameters by mapping the mean field deformation energy surface on the bosonic 
one is discussed and the applications listed. A critical assessment of the models is given 
with focus on the limitations due to the adiabatic approximation. The Tidal Wave ap¬ 
proach and the Triaxial Projected Shell Model are presented as practical approaches to 
calculate spectral properties outside the adiabatic region. 

Keywords: microscopic Bohr Hamiltonian; quadrupole excitations; triaxiality; IBM-1; 
triaxial projected shell model, tidal wave approach. 

21.10.Re, 21.60.Ev, 21.60.Fw, 23.20.Lv, 24.10.Cn 


1. Introduction 


Nuclei with mass number A > 50 are compact systems with a well defined sur¬ 
face, the shape of which and its orientation in space constitute collective degrees 
of freedom, which are discussed in any textbook on nuclear structure. This article 
reviews the descriptions of the quadruple mode and tries to asses, to which extend 
present nuclear theory is capable of describing and predicting the properties of the 
corresponding quantum states on the basis of the underpinning degrees of freedom 
of the nucleonic constituents. There are two major phenomenological descriptions 
of the collective quadrupole modes at low spin: the Bohr - Hamiltonian (BH), and 
the Interacting Boson Model (IBM), which describes the quadrupole mode in terms 
of boson operators that form a closed SU(6) Lie algebra. Both approaches are well 
exposed in textbooks, as for example.^®^ Section 2.1 reviews some basics of the 
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description in terms of the BH and recent developments of this phenomenology, 
where the aspects triaxiality are exposed in some detail. A microscopic version of 
the BH has been derived in the framework of all mean field approaches used in 
practice at present. Section gives an overview of derivations in the framework of 
the various Adiabatic Time-Dependent Mean Field approaches. Section covers 
corresponding work in the framework of the Generator Coordinate Method. Ap¬ 
plications to concrete nuclei are listed and selected examples discussed. Section 
presents some basics of the IBM phenomenology and the recent successful calcu¬ 
lation of its parameters by means of mapping the potential energy surface of the 
mean field. 

Nuclei are composed of a relative small number number of nucleons compared 
to other many-body systems. As a consequence, the ’’granular structure” of the col¬ 
lective degrees of freedom appears already after the excitation of few quanta, which 
results in a progressive decoherence of the collective modes. The left side of Fig. 
illustrates the point in a schematic way for a vibrational nucleus: The multi-phonon 
excitations encounter very soon the region of the quasiparticle excitations to which 
they couple. The right side shows a realistic Shell Model calculation. In accordance 
with the white adiabatic region on the left side, the collectively enhanced transi¬ 
tions are restricted to the one and two phonon states and higher up to the yrast 
region. The purely collective models covered in sections [^-[^ assume adiabaticity of 
the collective motion explicitly or implicitly. Their realm is restricted to the white 
adiabatic area in Fig. Outside, the coupling between the quasiparticle and col¬ 
lective degrees of freedom must be taken into account in a non-perturbative way. 
Section [^discuss the Tidal Wave approach and the Triaxial Projected Shell Model, 
which take this coupling into account. Section tries to assess the presently used 
models and points out some challenges, which seem important from the author’s 
point of view. 


2. The Bohr Hamiltonian 


The Bohr coordinate system and Hamiltonian are reviewed in detail by Prochniak 
and Rohozinski.l^ The collective states are represented by wave functions of the 
components (/i = —2, ..., 2) of the scale-free quadrupole deformation tensor 
which are expressed by the five-dimensional spherical polar coordinates 


a^ = j3 


>( 2 ) . 1 




.( 2 ) 


where ft are the Euler angles specifying the orientation of the shape. 


( 1 ) 


2.1. The Geometric Collective Model 

The Geometric Collective Model (GCM) is a parametrized version of the BH based 
on an expansion into scalars of increasing power constructed from the coordinate 
and their conjugate momenta tt^, which was introduced by Gneuss and Greiner .1^ 
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Fig. 1. Left: Schematic representation of the location of the collective quadrupole vibrational 
excitations relative to the quasiparticle excitations. The darker shades approximately indicate 
higher densities of quasiparticle states. From Ref.^ Right: Shell Model Calculation for ®^Ni. The 
width of the bars is proportional to the B{E2) value of the connecting transition. From Ref.^^ 


Only the quadratic term in tt^ is kept. In terms of the quadrupole deformation 
variables /3 and 7 and Euler angles O, the Bohr Hamiltonian is given by 


where 


and 
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The operator appearing in brackets in the kinetic energy is the Laplacian in five 
dimensions. Its angular part is the Casimir operator for the five-dimensional 
rotation group SO(5), which contains the rotations in physical space, acting on 
the Euler angle coordinates, as an SO(3) subgroup. The potential energy 
must be periodic in 7, with period 120 °, and it must be symmetric about 7 = 0 ° 
and 7 = 60°. Prochniak and RohozinskP discuss in detail the various methods for 
calculating the eigenvalues and eigenfunctions of the BH. 

The properties of the experimental 2^, 4^, 6 ^, 2^, 3^, 4^, Oj, 4^ states can 
be classified by assuming that the collective potential contains only three terms. 


Va{P, 7) = -^C2p^ - cos37 + (5) 

This provides a more complete scheme than the traditional classification into ro¬ 
tational and vibrational nuclei. CapricP^l demonstrated that the structure of the 






























June 23, 2015 0:26 WSPC/INSTRUCTION FILE ws-ijmpe-Frauendorf 


4 S. Frauendorf 








Fig. 2. Illustration of the qualitatively different shapes of the GCM potential function § obtained 
for different ranges of values for the parameter d. Potentials are shown as a function of /3, where 
7 = 0 (^ < 0 is equivalent with /3 > 0, 7 = tt/S). For / in MeV, the energy scale is also in MeV. 
Reproduced from Ref.E^ 


collective wave function is determined by two dimensionless parameter ratios only, 
the control parameter d = 112 ( 72 ( 74 /( 9 ^/ 5 ( 73 ) of the potential and the structure 
parameter h^S = T?C\I which controls the zero-point energy. The scale of 
the deformation parameter (3 is fixed by the ratio e = The scale of the total 

energy is ^ {\/be^B 2 ). For illustration, it is useful to introduce the scale of the 

potential energy / = = C^jC^. Then the structure parameter S = f /g is the 

ratio between the scales of the potential and the total energy. It controls to which 
extend the zero-point fluctuations of (3 and 7 wash out the details of the potential. 
In scaled form, the GCM Hamiltonian reads 


Hgcm2 = 9 


+ ^ + 


75 

h^S 




13'^ cos 37 + -p 



(6) 


where ^ = (3/e is the scale-free deformation parameter. The matrix elements of the 
charge quadrupole moments, which generate the E2 7 -transitions and the static elec¬ 
tric quadrupole matrix elements, are modeled by a homogeneously charged droplet. 




Jtt V ^ 



[a x a 


( 2 )* 


( 7 ) 


which fixes the scale e. 

Fig-i illustrates the different types of potentials. For 0 < d < 1 the potential 
has two minima at prolate deformation. For d < 0 there is a saddle at oblate shape, 
which connects smoothly with the prolate minimum via the 7 degree of freedom. For 
d —oc the limit of 7 -independence is approached. Caprio carried out a qualitative 
analysis of the ground state wave function by means of the WKBA. Comparing the 
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Fig. 3. Map of the GCM (d, S) parameter space. The regions in which qualitatively different 
structures occur are indicated. The curves S = S'<, and S = *S'> (dotted lines) provide estimates 
for the approximate boundaries between these regions. The squares show the location of ^°^Pd and 
^^^Sm in the parameter space. Bars along the edges of the plot represent structures which occur 
in their ideal form at d —)■ 00 or at S' = 0 or 00 . ”W-J” denotes the Wilets-Jean rigidly-deformed 
7 -soft structure. Reproduced and modified from Ref.^ 

zero-point energy with the extrema of the potential leads to a classification of 
structure of the low-lying collective states, which is mapped in Fig. 

At the line 5'>, the zero-point energy is equal to the height of the barrier between 
the two axial minima. Above this line the wave function spreads over both minima. 
Below the line 5'>, the wave function becomes progressively suppressed under the 
barrier. That is, the region of deformed nuclei lies within the rectangle [d < 1, 5'>]. 
At the line 5'<, the zero-point energy is equal to the energy difference between the 
oblate saddle and prolate minimum. That is, the prolate nuclei are located below 
the line 5'<, which is denoted by ”rotor-vibrator” region. The ’’stiffness” for the {3- 
and 7 -vibrations varies with d (double arrow), where ’’stiffness” refers to the relative 
order of the Oj and 2^ states. The lines 5'> and 5'< demarcate the region of 7 -soft 
nuclei, for which the wave function extends over the whole range 0 < 7 < 7 r/ 3 . 
Outside the rectangle [d < 1 , 5'>], the wave function changes gradually from a 
harmonic oscillator to a quartic oscillator, for which the term confines the wave 
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function. lachellof^ introduced two schematic potentials that classify the structures 
in this region. The ”E(5)” potential does not depend on P and 7 for < Pw where 
it jumps to 00 . The wave functions of the lowest states resemble the ones in the 
region of large values of —d < —5, where the potential weakly depends on 7 . The 
”X(5)” potential has a term oc 7 ^ added to the E(5) potential. The corresponding 
wave functions resemble the ones in the region — 5 < d < 1, where the potential 
prefers the prolate shape. It has become popular to use the acronyms E(5) and 
X(5) for these structures. 





Fig. 4. Experimental level scheme and B(E2) strengths for The number in parenthesis 

under the transition energies (keV) are the B(E2) values (e^ fm"^) for the transitions. Data from.^J 
Preparation of the figure by A.D. Ayangeakaa is acknowledged. 


Eig.j^shows ^^^Pd as an example for an E(5) nucleus .^Eitting the experimental 
energies and E 2 transition probabilities provides the GCM parameters d = —43 and 
= 24 X 10“^^ = 56 X 10^^(MeVs)“^d^, which corresponds to 5'/5'> = 0.9. The 
GCM predictions are shown in Eig.|^ The agreement with the data is characteristic 
for GCM phenomenology. The GCM predictions are similar to the ones of the E(5) 
model (e. g. Eig. 8 of Ref.^. An example for the X(5) structure is ^^^Sm.C^The 
GCM parameters are d = —1.75 and fi^S = 71 x 10“^^ = 18 x 10^^(MeVs)“^d^, 
which corresponds to S/Sy = 0.25. The ratios of the energies and B{E2) values 
from the X(5) model and from the GCM turn out to be similar as well (c. f. RefP). 
The reason of the similarity of the GCM fits and the schematic models is that the 
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Fig. 5. GCM predictionJi^ of the energies and B{E2) strengths for ^®^Pd normalized to the 
experimental E(2'^) and and B(E2]2^ —)► 0^) values. The number in parenthesis under the 
transition energies (keV) are the B{E2) values (e^ fm"^) for the transitions. Preparation of the 
figure by W. Li is acknowledged. 


wave functions of the considered few lowest states do not resolve the details of the 
potential, because the wavelength is larger. CapriJ^ presents figures of the ratios 
E{I)/E{2^) and B{E2;I E)!B{E2;2^ ^0^), which allows one to extract the 
GCM parameters, and discusses fitting strategies. A GCM analysis of the = 90 
region can be found in his thesis,ini where he also analyses the difference between 
the GCM and the schematic models. 


2.2. Triaxiality 

Strong deviations from axial shape have been in the focus of recent research. The 
two-parameter version of the GCM encompasses the range between the axial regime 
with a stiff y-vibration and Wilet-Jean (WJ) limit of a y-independent potential. To 
account for a stabilization of the triaxial shape, a term with a minimum at y = tt/G 
has been added to the potential ([^, where the authors used different functions 

CapricP^ studied the consequences of the gradual stabilization of the triaxial 
shape in a systematic way for a frozen deformation /3 = Po hy diagonalizing the 
Hamiltonian 


HgCM3 = + X [l - cos 37 + ^ cos^ 37 ] . 


(8) 
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As already noticed by Zamfir and Castenf^ the order of the two signatures of the 
quasi-gamma band ( 2 ^, Sf, 4^, 5^, 6 ^, ..) provides a clear spectral significance for 
the type of triaxiality. The order is easily understood considering the two limiting 
cases. The y-independent potential (WJ) has SO(5) symmetry. The states 2 ^, (3^, 
4^), (5]^, 6 ^), ... belong to multiplets with SO(5) seniority v=l, 2 , 3, ... , where the 
energy is v(v+3)/4 in units of the ^(2^). That is, the even-/ states are below the 
average of their odd-1 neighbors. When the triaxiality parameter is fixed at 7 = tt/G 
(TR), the ratio of the three moments of inertia is 1/1/4 (irrotational flow). The 
energies for this triaxial rotor (TR) are given by /(/ + 1)/12+ [21 {n + 1/2) — n^] /4 
in units of the E{2f)^ where n is the number of wobbling excitation and I is even 
(odd) when n is even (odd). That is, the odd-/ states are below the average of their 
even-/ neighbors. In the case of an axial rotor (AR), the two signatures of the 7 
band are degenerate. 

The staggering parameter 

S{I) = {E{I) - 2E{I - 1) + E{I - 2 ))/E( 2 +), (9) 

calculated from the energies E{I) of the quasi-gamma band, is used to characterize 
the triaxiality .113 The limiting models give for, respectively / even/odd, 

Swjil) ^ 8 T 4 + 2 ) ’ ^ 6 ^ “ 2 ) ’ ^ 3 ’ 

Another spectral significance is the ratio /^(2j)/E'(4^) which is > 1 , 1 , < 1 for the 
AR, WJ, TR, limits respectively. CapricP^ studied in detail the transition between 
the regimes in the framework of the model <§ and provided the figures and tables of 
the ratios of the energies and of B{E2) values that allows one to determine the model 
parameters. Fig. [^illustrates in the upper part the development from the AR to the 
WJ regime and in the lower part the development from the WJ to the TR regime. 
It is noted that small staggering is not necessarily an indication for the common 
case of small triaxiality. It also appears when the staggering phase changes sign at 
the transition from the WJ to the TR regime. The low ratio /^(2j)//^(4^) < 1 can 
be used to discriminate this case from axiality. 

As an example. Fig. [^demonstrates the change of triaxiality along the Pd iso¬ 
tope chain. The ratio /^(2j)//^(4^) stays around 0.85, with a very shallow dip 
around N= 68 , which indicates strong deviations from axiality. Below N = 66 one 
sees the y-soft pattern: even I below odd /. For N = 68 and 70, the triaxial rotor 
pattern emerges for S{I > 5): odd I below even /. There is a vague indication that 
the y-soft pattern returns for N = 72. The Ru chain shows a similar trend as the 
Pd chain, where the amplitude of the staggering is larger: S{8) — S(9) = 1.2 in 
^^^Ru and S{8) — S(9) = 0.6 in Pd. McCutchan et al^ studied the systematics 
of the staggering parameter. Fig. [^ from their work collects nuclei with the odd-I- 
low staggering of a stabilized triaxial shape. The other investigated nuclides either 
show the ^'-positive pattern of axial shape or the even-/-low pattern of y-soft nuclei. 
Toh et reported odd-/-low staggering for ^^Ge, where the neighboring nuclides 
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Fig. 6 . Staggering of level energies within the quasi -7 band, as measured by the staggering pa¬ 
rameter S{L) calculated for the Hamiltonian § with = 0 , for (d) x = 50, (e) x = IOO 5 cind 
(f) X = 200. in the upper part and for (d) ^ = 0 with x = 50, (e) ^ = 0.5 with x = 100, and (f) 
^ = 0.8 with X = 500 in the lower part. The potential V('y) is shown in the inset, with the ground, 
quasi- 7 , and quasi -77 band head energies indicated. Reproduced and adapted from Ref.^l^ 


show the even-/-low pattern. However this cannot be taken as ’’Evidence for rigid 
triaxial deformation at low energy” as claimed by the authors, because the observed 
value of S{S) — S{9) = 0.5 is an order of magnitude smaller than the rigid rotor 
value of 12 given by Eq. (10). The nucleus seems to be located in-between the cases 
e) and f) shown in the lowest panel of Eig. That is, the zero-point energy is sub¬ 
stantially larger than the depth of the potential minimum at triaxial shape, which 
is far from a rigid-triaxial regime. The authors carried out Shell Model calculations, 
which well reproduce both experimental energies and B{E2) values (amplitude of 
the staggering in particular). This is expected because the Shell Model takes the 
dynamics of the quadrupole degree of freedom into account. 
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Fig. 7. Experimental staggering parameter S{J) for the Pd isotopes. The black line shows the 
ratio F;(2+)/F;(4+). Data from ENSDhEHand RehEHl 


3. Adiabatic Time-Dependent Mean Field approach 

The Adiabatic Time-Dependent Mean Field (ATDMF) approach is the standard 
method for deriving the parameters of the BH from the Fermionic underpinning. 
Baranger and Kumar introduced it in their pioneering work.l^ ADTMF is very well 
exposed in the reviews by Prochniak and RohozinskP and Niksic, Vretenar and 
RingHH Here we only sketch few steps that are essential for the discussion to be 
followed. 

First, the potential is obtained by the mean field approach of choice, which 
is constrained to provide a given value of the expectation value of the microscopic 
mass quadruple moments go = (Qo) and g 2 = (Q 2 )- The deformation parameters go 
and g 2 and the three Euler angles, which specify the orientation of the quadrupole 
shape, are the dynamic variables of the BH. To keep contact with phenomenology 
one often introduces the standard dimension-less deformation parameters P and 7 
by the relation 

-A- (^roA^/^j , ro = 1.2/m, (11) 

which assumes the liquid drop relation between the deformation parameters and 
the quadrupole moments. Second, a classical BH is derived by time dependent 
perturbation theory. The adiabatic approximation is used, which means that only 
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minima at odd J 



Fig. 8. Experimental staggering parameter^ J) for nuclides that shows the staggering associated 
with triaxial shapes. Reproduced from Ref.^^ 


terms quadratic in the time derivative of the expectation values (Q^) are kept. 


Hclass — '^vih T ^rot T T)? 

(12) 

Tvib = 

(13) 

^ Ji = 4.Bip‘^ sin ^(7 - ^7ri), 

(14) 


i=l,2,3 


where uji are the angular velocities with respect to the three body fixed axes. The 
classical BH corresponding to the GCM Hamiltonian <§ is obtained by setting 
= Bi = ^/b/2B2 and B^^ = 0. The microscopically calculated mass 
parameters are far from this drastic simplification (see below), which is a severe 
limitation of the GCM phenomenology. 

Third, the classical BH (12) is quantized. The procedure is involved because 
the mass parameters depend on the deformation. It is described in the reviews^^ 


where the resulting complicated expression for the quantal BH is quoted and the 
numerical methods for the solution of the eigenvalue problem are discussed. The 
calculation of the potential and the mass parameters depends on the mean field 
approach of choice. It has recently become customary to refer to the various variants 
as 5DBH (5 Dimensional Bohr Hamiltonian). 
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3.1. Pairing plus Quadrupole Model 

The generic version of the microscopic 5DBH-PQQ, is derived from the 

Pairing plus Quadrupole-Quadrupole (PQQ) Hamiltonian, 

Hpqq = - GP^P - 1 ^ gt (15) 

which combines a spherical term hsph constructed from adjusted spherical single 
particle energies with a Monopole Pairing interaction for the short-range correla¬ 
tions and a Quadrupole-Quadrupole interaction for the long-range correlations. It 
generates the mean field Hamiltonian 

hpQQ — hsph T A P^ P fqoQq + hvq2 {Q 2 + Q— 2 ) ~ XN. (f^) 


This is the BCS Hamiltonian for a deformed Nilsson-type potential, which has the 
standard BCS ground state |) = | A, A, go 5 ^ 2 )- (Protons and neutrons are not explic¬ 
itly distinguished for simplicity.). The term —XN assures the correct expectation 
value of the particle number by the equation {N) = N. The pair field is deter¬ 
mined by the selfconsistency equation A = G (P). The pair field A (go, Q 2 ) and the 
chemical potential A(go,g 2 ) are functions of the deformation parameters. 

The mass parameters take the Inglis-Belyayev (IB) form with /i = 0, 2 


M 


(n) ,fiiy 


= E 


kl 


= 


{Hh 


-M 




Ji = 


{\Q^kl) {kl\Q^\) 

{Ek+EiT 

P\Ji\kl) {kl\Ji\) 


kl 


{Ek + El) 


(17) 


(18) 


where |) and \kl) stand for the zero- and two-quasiparticle states, and Ei are the 
quasiparticle energies. The term P^^, takes the deformation dependence of A and 
A into account. The somewhat complicated expression is given in Ref.l^ The mass 
parameters for the dimension-less deformation variables P and 7 are obtained by 
transforming the matrix as cP( 7 )^PcP( 7 ), where cTl{j) is a rotation by 7 
and a multiplication by the scale factor c. 

The PQQ model takes two shells into account. The pairing strength is adjusted 
to the experimental even-odd mass differences. The quadrupole coupling constant 
p is treated as a parameter that is adjusted to reproduce the experimental spectra. 
The IB mass parameters turn out to be to small. A common scaling factorP^ ^ 2 to 
3 is multiplied to obtain the experimental energy scale. For the E2 matrix elements 
a polarization charge of kZjA with k ^ 1.6 is used. The two parameters set the 
scales for the energy and charge quadrupole moments. Compared to the GCM 
and IBM phenomenology, which use the same free scales, the 5DBH-PPQ model 
has only one parameter, k, which changes smoothly (x through regions of 

major re-structuring. The nature of collective excitations is determined by the shell 
structure of the nucleonic orbitals. In case of the PQQ, it is encoded in the spherical 
single particle levels, which are input as well, though determined from data different 
from the collective excitations. 
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Fig. 9. The potential V (left), the mass coefficient (middle, arbitrary units) and the ground 
state wave function (right) calculated by means of the 5DBH-PPQ for ^^^Xe. Reproduced with 
permission from Ref.^^ 


Kumar, Baranger, and Gupta applied the 5DBH-PQQ, which they call Dy¬ 
namic Pairing Plus Quadruple (DPPQ) model, to nuclei with (Z,A)= (74,182- 
186), (76,186-192), (78,192) in Ref. ,>221 (72,166-168) in Ref.,l231 (62,150-152) in Rei.^ 
(62,142-154) in ReL^ (62,150) in Ref.l22l (64,154-156) in Ref.,l2Hl (66,154-156) 
in Ref.,l2^ (60,144-150) in Ref.,1^ (58,130-136) in Ref.,ISl (56,122-134) in Rei.^ 
(54,124) in Ref. ,1^ and the N = 66 isotones with Z = 42 — 58 in Ref.l^In addition 
to the energies and E2 matrix elements, they calculated magnetic moments, spec¬ 
troscopic factors, and EO moments. They investigated the phase transition region 
around A/" = 90 in considerable detail, subsequently tuning the model parameters. 
More recent work extensively compares the 5DBH-PQQ with the IBM phenomenol¬ 
ogy* 

As an example, consider ^^^Xe. 5DBH-PQQ well reproduces both energies and 
B{E2) values. Eig. illustrates the calculations. The potential has a minimum at 
axial shape around ;d=0.24. The mass coefficient 5^/3 is far from being constant, as 
assumed in the GCM phenomenology. It increases toward small /3 on the prolate 
side and becomes large on the oblate side. The wave function tends to concentrate in 
regions of large mass like in regions of low potential. This has the consequence that 
the wave function has its maximum around /3 = 0.16 and does not change much in 
7 -direction. That is, it has the character of the WJ model. The calculated staggering 
of the quasi -7 band is of the even-/-low type in agreement with the experiment, al¬ 
though the model overestimates the experimental staggering parameter: S'(5) = 208 
keV and S{6)= -92 keV to be compared with the calculated values 498, -475 keV, 
respectively. The ratio B{E2,3^ 2+)/5(E2,3+ ^ 2+) is small (0.024 experi¬ 

mental and 0.064 calculated) whereas the ratio 5(E2,3+ ^ 4+)/5(E2,3+ ^ 2+) 
is large (0.26 experimental and 0.19 calculated), which reflects the seniority selec¬ 
tion rules of the WJ limit corresponding to 0 and 0.4, respectively. The calculated 
energy of the 7 band head E(+) = 1097 keV is larger than the experimental value of 
847 keV and the calculated E(0j) = 1099 keV is lower than the experimental value 
of 1269 keV. The example has the typical accuracy the 5DBH-PQQ calculations. 
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Fig. 10. The experimental level scheme of (black levels) compared with the ones calculated 

by means of the 5DBH-SC model (red levels). Reproduced from Ref.^^ 


It demonstrates that the just looking on the calculated potentials for guessing the 
structure of the low-lying collective excitations may be risky, because the deforma¬ 
tion dependence of the mass coefficients may cause substantial modifications of the 
wave functions. Conversely, a potential derived by fitting data in the framework of 
the GCM, which assumes a simple deformation dependence irrotational flow for the 
mass coefficients, may substantially differ from the microscopic potential calculated 
by the mean field approach. 


3.2. Pairing plus Quadrupole Model with Local Random Phase 
Approximation 

Matsuo, Matsuyangi, Nakatsukasa et a/EEll removed the problem with the too 
dilute energy scale of the 5DBH-PQQ approach by adding a quadrupole pairing 
interaction to the PQQ Hamiltonian (15). Quadrupole pairing is known to increase 
the rotational moment of inertia by about 30% when calculated by means of the self- 
consistent cranking model. The pair field generated by the quadrupole pairing adds 
a new term to the IB value, which increases it to the Thouless-Valatin (TV) value. 
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Fig. 11. The experimental invariants (Q^) = 5 ^ and cos 87 = cos 87 ) / for 

compared with the ones calculated by means of the 5DBH-SC model (QCBH dyn). Reproduced 
from Ref.^ 


The Adiabatic Self-consistent Collective Coordinate (ASCC) method developed by 
the authors allows them to calculate the TV correction to the IB values for the mass 
parameters and as well. In practice, the equations of the Quasiparticle 
Random Phase Approximation (QRPA) are solved for each point of the grid 
(Local QRPA -LQRPA). The mass parameters are derived from the frequency and 
oscillator length of the lowest two solutions in the standard way. An increase of 
20% to 30% of all mass parameters of the BH is found, which results in the energy 
scale of the experiment. 

Carrying out the LQRPA for the large number of grid points is computational 
extensive. For this reason only nuclei with mass below 100 have been studied: 
(Z, N) = (34, 34 - 38) in Ref.,1^ (24, 34-44) in Ref. pHEH and (12, 18-22) in Ref.SSl 
Without any free parameters, the agreement of the calculations with the experiment 
is comparable with the 5DBH-PQQ, which uses a scaling factor. The analysis in 

Refimnol 

demonstrates that the two-parameter BH phenomenology in Sect. |2.l| is 


not capable of adequately accounting for the microscopy based wave functions. In 
the case of Oj states, the deformation dependence of the mass parameters leads to 
structures in between the /3-vibration and shape coexistence. 
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3.3. The Micro-Macro Model 

Rohozinski et introduced the 5DBH-SC version of the microscopic BH. The 
potential energy V (/3, 7 ) is obtained by means of Strutinski’s Shell Correction (SC) 
method, which combines the deformation energy of a droplet of nuclear matter 
with a shell correction obtained from the single particle levels in a deformed mod¬ 
ified oscillator (Nilsson) potential. The SC method, which is now often called the 
Micro-Macro approach, is well tested for calculating deformation energies. It uses 
the spherical single particle energies and the surface and Coulomb energy of the 
nuclear droplet as input, which are determined from data other then the low-energy 
collective excitations. The scale of the E 2 matrix elements is fixed by the Coulomb 
field of the droplet. The mass parameters are calculated by means of the IB ex¬ 
pressions (18), which are too small when the standard BCS monopole pairing is 
used. Prochniak et al^ removed the problem by improving the treatment of the 
pair correlations. They construct a collective Hamiltonian for the pair gap A and 
the gauge angle ^ in analogy to the BH for the quadrupole degrees of freedom. For 
each point on the P-j grid, they calculate the pairing ground state wave function 
and determine the most probable value of A, which they the use in calculating the 
IB mass parameters. The most probable value of A is smaller than the BCS value 
and depends less on deformation, which brings the mass parameters to the right 
scale. With this modification the 5DBH-SC becomes parameter-free. ProchnialP^ 
demonstrated that including the direct coupling between the quadrupole and the 
pairing degrees of freedom in the framework of a nine dimensional generalized BH 
provides essentially the same results as the approximatioiP^ of using the the most 
probable A values. The fact that the right scale of the mass parameters can be 
achieved either by the dynamical treatment of the monopole pairing or by inclu¬ 
sion of the quadruple pairing (see preceding Sec. 3.2) raises the question about the 
consequences of combining both improvements of the BCS treatment. 

The 5DBH+SC has been applied to nuclides with Z, N) = (52-62, 68-80) in 
Ref.^ (46, 60-64), (44, 60-70) in Ref.,l^and (44, 60) in Ref.®^ As an example. Fig. 


10 demonstrates that the 5DBH+SC very well reproduces the spectrum of ^^^Ru 


without fitting any parameters. The 5DBH+SC has been preferentially applied to 
nuclei for which the E2 matrix elements have been measured by means of Coulomb 
excitation, providing a stringent test, which is well met by the model. In addition 
to direct comparison, the authors studied invariants of the quadrupole mode, which 
are derived by combining measured transition matrix elements. Srebrny and Cline 
describe the method in Ref. where they discuss a number of examples and expose 
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the history. It exploits the sum rules 

o2 ^ 

^ (/Ja|/32|/Ja) = ( 2 ^ 21 ^ 100 ) (/Ja|Q^Q,|/Ja) 

= Y, {2M2i^|00) (19) 

gul' M' a' 

jj^\ _ 1^pl2v\2k) ,iJla\Qy,QyQ\\lla) 

gvnX 

= ^ {2ii2u\2n) (2 /i:2 A|00) 

livKXI'M'a'I"M"a" 

X {IIa\Q^\rM'a') {!'M'a'\Q,\I"M"a") {!"M"a"\Qx\IIa ), (20) 

which are evaluated with the experimental matrix elements. The theoretical ex¬ 
pectation values are directly calculated by integrating the wave functions. Fig. 
exemplifies the accuracy of the model for ^^^Ru. The calculations overestimate the 
increase of the deformation with I within of the ground band. The nucleus is 7 -soft 
with a preference of prolate shape, which is consistent with its even-/-low staggering 
of the 7 band. 

3.4. Mean Fields on a Microscopic Basis 

More recent versions of the microscopic BH are based on mean fields derived from 
the Skyrme Energy Density Functional (SK-EDF), the Relativistic Mean Field ap¬ 
proach (RMF), and the Hartree-Fock-Bogolybov method applied to the Gogny ef¬ 
fective interaction (GI). Gompared to the PQQ and SC versions of 5DBH, the 
spherical single particle energies are no longer phenomenological input (either di¬ 
rect or as the parameters of the Nilsson potential), but derived from a more funda¬ 
mental layer. However this does not necessarily imply that they reproduce the shell 
structure more accurately than the SC or PQQ approaches. The detailed Z—, N- 
dependence of the binding energies is excellently reproduced by the Finite Range 
Droplet Model (FRDM) by Moller and collaborators, who base the SC approach on 
the folded Yukawa potential (c. f. Ref.E^and references therein). The spectral prop¬ 
erties obtained by the 5DBH are very sensitive to the shell structure, and therefore 
to the accuracy the mean field reproduces it. In view of this, an application of the 
5DBH-SC method to the FRDM promises a relatively simple and computational 
inexpensive way of predicting the collective quadrupole excitations of even-even 
nuclei. 

Deriving the 5DBH from the modern mean field approaches requires several 
sophistications, which are well exposed in the review articlesElEll The definition of 
the collective coordinates is less direct than for the PQQ and SC versions. The 
collective coordinates ca^{\ 2 y) = {Tnf^\ 2 v\Q^\mf^\ 2 v) are implicit functions of 
the Lagrange parameters of the constraints —\ 2 yQv that are used to generate the 
manifold of mean field states |m/, \ 2 y)‘ This leads to a modification of the IB mass 
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^M\ 


«), 


The interaction induces 


parameters of the form 

a time-odd terms in the mean field, which generate corrections to the IB values 
of the mass parameters increasing them to TV values. Since their evaluations is 
difficult (see Sec. 3.2), they are usually neglected or only partially included. Further, 
there are zero point energy corrections (ZPE), which can be derived in different 
ways. Usually the expressions obtained from the Generator Coordinate Method with 
the Gaussian Overlap approximation (GCM+GOA, see below) are used. The ZPE 
consist a rotational part AVrot = -i i and a vibrational 


part AVvib = Tr (a1(2)AJ(3)) /4. 

Prochniak et alW^ introduced the 5DBH+SK, which is based on the Skyrme 
EDE. They use the IB mass parameters, no ZPE, and either a monopole pair 
interaction or a zero-range pair interaction. The mass parameters turn out to be 
too small, which is corrected by scaling the energies by a factor of 1.2 -1.3. They 
studied the nuclides with (Z, N)= (40-42,62), (44,66), (56, 64), (54-56,70) inRef.,1^ 
(92-94,146), (96,150-152), 98,152-154) in Ref. ,>^3 (48, 52-58) in Ref.,Eni (42, 42-58) 
in Ref.,ISIl (42, 54-58) in Rel^ (42, 58) in Ref.,Eniand (36, 36-40) in Ref.ESiLike the 
5DBH+SC, the 5DBH+SK has been preferentially applied to nuclei for which the 
E2 matrix elements have been measured by means of Coulomb excitation. 

The Bejing-Munich-Zagreb collaboration worked out 5DBH+RME, which is 
based on the Relativistic Mean Eield approach. It is well exposed in the review .1^ 
It uses IB mass coefficients and ZPE corrections. The pairing is taken in mean field 
approximation for a density-dependent (^-interaction or a separable in momentum 
space interaction. In most cast cases, the IB mass coefficients give a too diluted 
spectrum, which is corrected by a scaling factor. An efficient method to calculate 
the TV contributions has been suggested in Ref. ,1^ however not yet implemented. 
They studied the nuclides with (Z, N)= (94, 156), (68, 98) in Ref.,1^ (60, 84-94), 
(62, 88-90), (64, 88-90) in Ref.,l^(60, 90) (sensitivity to pair correlations) in Ref.,1^ 
(54-56, 74-80), in Ref.,1^ (50, 52-80) in Ref. ,1^21 (38-40, 60) in Ref. ,1^11 (36, 32-50) in 
Ref.,1^ (30, 34-38) in Ref.,l^and (14-18, 28) in Ref.^^ 

Delaroche and collaborator^^ used the 5DBH+GI, which is derived from the 
Gogny Interaction. By means of self-consistent cranking for rotation, they include 
the TV corrections into the moments of inertia. The vibrational mass coefficients 
are of the IB form. ZPE corrections are taken into account. There is no scaling of 
the energies. In a bench mark study, they carried out calculations for all even-even 
with 10 < Z < 100 and 200 < 200. The results for the energies and E2 and EO 
matrix elements for the yrast levels with / < 6, the lowest excited 0+ states, and the 
two next yrare 2+ states accessible in the form of a table as supplemental material 
to the publication. A thorough statistical analysis of the merits of performance has 
been carried out. The authors state: ”We assess its accuracy by comparison with 
experiments on all applicable nuclei where the systematic tabulations of the data are 
available. We find that the predicted radii have an accuracy of 0.6%, much better 
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than can be achieved with a smooth phenomenological description. The correlation 
energy obtained from the collective Hamiltonian gives a significant improvement to 
the accuracy of the two-particle separation energies and to their differences, the two- 
particle gaps. Many of the properties depend strongly on the intrinsic deformation 
and we find that the theory is especially reliable for strongly deformed nuclei. The 
distribution of values of the collective structure indicator RA2 = ^(4^ )/E{2 ^) has a 
very sharp peak at the value 10/3, in agreement with the existing data. On average, 
the predicted excitation energy and transition strength of the first 2 + excitation 
are 12 % and 22 % higher than experiment, respectively, with variances of the order 
of 40-50%. The theory gives a good qualitative account of the range of variation of 
the excitation energy of the first excited 0 + state, but the predicted energies are 
systematically 50% high. The calculated yrare 2+ states show a clear separation 
between 7 and /3 excitations, and the energies of the 2 + 7 vibrations accord well 
with experiment. The character of the 0 J state is interpreted as shape coexistence or 
/3-vibrational excitations on the basis of relative quadrupole transition strengths. 
Bands are predicted with the properties of (3 vibrations for many nuclei having 
RA2 values corresponding to axial rotors, but the shape coexistence phenomenon 
is more prevalent.” In addition they observe that the Oj states are generally ”too 
vibrational”. 


4. Generator Coordinate Method 


The Generator Coordinate Method (GCM - unfortunate coincidence with GCM 
for the Geometric Collective Model) constructs the collective wave function as a 
superposition of the constraint mean field solutions |T) = J d^af{a)\mf, o), where 
a = {o_ 2 , ..., 02 } is a short hand notation for the five components. Minimizing the 
energy leads to the Hill-Wheeler (HW) eigenvalue problem for the weight function 

/(a) 

j d^a' {a\H\a’) f{a') =E j d^a' {a\\a’) f{a'). ( 21 ) 


The overlap kernel ((a| in the integral equations appears because the set of mean 
field solutions ja) represents a non-orthogonal basis. 

The Gaussian Overlap Approximation (GOA) by Girod and GrammaticoJmUll 
allows one to recast the integral equations into the form of the standard 
5DBH differential equation. It approximates the overlap kernel by a Gaussian 


( 0 :^ 117 ) = exp 
energy kernel {a 


Tij.u9{dt)iiv {(Xf, - 7 ) > where a = {a + a')/2, and the 

H\a') by a second order Taylor expansion around {a\\a'). The 
GOA provides expressions for the ZPE. The mass coefficients are of the Peiers- 
Yoccoz (PY) type, which are known to be smaller than the IB values. For this 
reason they are replaced by the IB expressions when the GOA is used for mapping 
the HW integral equations on the 5DBH differential equations. 

The direct solution of the HW integral equations avoids the second order Taylor 
expansion of the energy kernel, which is one part of the adiabatic approximation. 
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Fig. 12. The spectrum of calculated by means of the 5DBH-GI.I^ The number in paren¬ 

thesis under the transition energies (keV) are the B{E2) values (e^ fm"^) for the transitions. 
Preparation of the figure by A.D. Ayangeakaa acknowledged. 


However using only the ground state configuration of the mean field solutions as 
HW basis assumes that the collective motion is slow enough that the coupling to 
excited mean field configurations can be neglected, which is another part of the adi¬ 
abatic approximation. The HW equations are reformulated using the deformation 
parameters [3 and 7 and the three Euler angles O. The integration over the angles 
takes the form of the projection operators on good angular momentum 
the integration of the deformation parameters a = {/d, 7 } is discretized. The HW 
equations become 

fxicij) = ( 22 ) 

jK jK 

Restriction to axial shapes 7 = 0 and M = K = 0 considerably reduces the 
numerical effort. Axial GCM calculations have become routine for the various mean 
field approaches, which cannot be reviewed for space limitation. The references can 
be found in the articles about the triaxial extension to be discussed below. An 
illustrative example is detailed study of ^^^Sm in the framework of the RMF in 

Refm 

The direct solution of the HW is numerically very demanding because of the 
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Fig. 13. Low-lying spectra and B{E2) values (in fm^) of Results from (b) the full 

RMF-GCM calculation compared with (c) 5DCH-RMF results and with (a) experimental data. 
Reproduced from 


dimension of five. For this reason it has been applied only to light nuclei so far. Ben¬ 
der and HeeneiP^ studied generating the HFB mean field solutions from the 

SK-EDF combined with a zero-range pairing interaction. The mean field solutions 
are projected on good particle number. The study focuses on the new aspects of 
triaxiality, like the mixing of the angular momentum components K with respect to 
the body-fixed frame. Rodriguez and EgidcP^ used particle number projected HEB 
solutions generated from the Gogny interaction to carry out calculations for ^^Mg. 
Comparing with the full triaxial calculation, they observe that the ground band and 
the P band on the 0 J state are rather well reproduced by the axial approximation. 

Yao et al. carried out a benchmark calculation for ^^Kr, which were based 
on the RME. Eig. p!3| displays the results, which are compared with 5DBH+RME 
calculations using the same interaction. The results agree rather well. The 5DBH 
spectrum is somewhat too diluted, while the GCM gives the right scale. This comes 
as a surprise, because earlier studies found that the PY moments of inertia of GCM 
are smaller than the IB cranking values. The reason is unclear. The low-lying states 
of ^^Kr have been interpreted in terms of a prolate ground state coexisting with 
an excited Oj state of smaller oblate deformation. The deformation difference is 
manifest by the spacing of the two rotational bands built on the band heads. Both 
approaches reproduce the experiment in considerable detail, which is also the case 
for the 5DBH+SK (Eig. 10 of Ref.l^ and 5DBH+GI (Eig. 16 of Ref.l^ calculations 
for this nucleus. Tab. II of Ref.l^ compares calculated spectroscopic quadrupole 
moments with experiment. The figures and the table represent the state of art of 
microscopic approaches based on self-consistent mean field approaches. 
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0 ( 6 ) 



U(5) X(5) SU(3) 


Fig. 14. The IBM parameters of selected isotope chains arranged in the Symmetry Triangle. 
The symmetry limits of the subgroups are indicated at the three corners. The ’’transition point 
symmetry” X(5) marks the transition from the vibrational regime (limit U(5)) to axial deformation 
(limit SU(3) and E(5) marks the transition to the 7 - soft deformed regime (limit 0(6)). Adapted 
from Ref.^^ 


5. The Interacting Boson Model 

The collective quadrupole mode is described by means of the creation operators for 
nucleon pairs with spins 0 and 2, which are called s- and d- bosons and denoted by 
and respectively P They form the closed Lie algebra of the SU(6) group. In 
this review only the simplest version, the ’’Extended Consistent Q Formulation” of 
IBM-1 of Warner and Casten,!^ ig presented. The Hamiltonian contains two IBM 
parameters and the energy scale, such that: 

HiBMAiC, X) = fl((l - Ond - ■ Q^), (23) 

where fid = and = [s ^The Hamiltonian is diagonalized 
within the space of fixed number of bosons, Nb = rig rid, which is taken to be 
half the number of valence nucleons. The matrix elements of the charge quadrupole 
moments are taken to be proportional to with an effective boson charge fixing 
the scale. As in the case of the GCM, the two parameters ( and y determine the 











June 23, 2015 0:26 WSPC/INSTRUCTION FILE ws-ijmpe-Frauendorf 


The Low-Energy Quadrupole Mode of Nuclei 23 


character of the collective states. 

It has become custom to map the structure of the collective states on the ” sym¬ 
metry triangle”, which is a polar plot of ( and y with the radius p and the angle 0 
defined by 

^ tt ^ ^ 271 


P = 


^/3 cos 0^ — sin 0^ ’ 


0 = 


0. 


X’ 


^x = ^ 

^ 3V7 


X- 


(24) 


The three corners of the triangle are parameter combinations that generate addi¬ 
tional symmetries with respect to the subgroups U(5), 0(6), SU(3), which corre¬ 
spond to the harmonic vibrator, the 7 -independent, axial rotor limits of the GCM 
approach. It is an attractive feature of the IBM that simple algebraic expressions 
describe the energies and reduced transition probabilities of the three symmetry 
limits. The two-parameter triangle is used to classify nuclei, where the concept 
of ’’Quantum Phase Transitions” is invoked. Order parameters Pb^b are defined 
by mapping the IBM Hamiltonian on a basis of coherent states which represent a 
condensate of d-bosons. 




where 


= 


5 + + 


/3b(cos(7b) 7 + + c^7))- 


V2 


(25) 


(26) 


The expectation value of the IBM Hamiltonian with the coherent state (25) has 
been given by Ginoccio and Kirson,!^ 

Eibm{Pb,1b) = {Nb, Pb,1b\Hibm\Nb, Pb,1b) 


= Ce 


/ 7c + ((1 - ONb - K(1 + x')) (/3 b) 

7(/Vb - 1)(/3b)7 

V l + {pB? 

V (1 + (/3b)2)2 ) 


(^1 - Y 7^xPbcos{3jb) + 


(27) 


Border lines between the regimes characterized by the symmetry limits are defined 
by applying the Ehrenfest classification of instabilities of thermodynamic poten¬ 
tials to the energy function E/bm(/^b,7b), which are shown in Fig. Only the 
thermodynamic limit of an infinite system leads to a sharp phase boundary. The 
instabilities of E/bm(/^b, 7 b) demarcate the center of a cross-over region between 
the three regimes. In order to describe the instabilities by simple algebraic expres¬ 
sion, lachellcP^ invoked the discussed solutions of the BH with a schematic square 


well potential (c. .f Sec. 2 . 1 ). These are commonly referred to as the X(5) and 
E(5) ’’transition point symmetries” although they are no algebraic symmetries that 
correspond to a symmetry group like the corners of the triangle. 

The two-parameter fits usually well account for the relative energies of the 2^, 
4^, 6 ^, 2^, 3^, 4^, Oj, 4^ states and the relative B{E2) for the transitions between 
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Fig. 15. IBM-1 fit to the spectrum of The number in parenthesis under the transition 

energies (keV) are the B{E2) values (e^ fm'^) for the transitions. The IBM parameters are ^ = 
0.61, X = 0.58, Nb = 5. Preparation of the figure by W. Li is acknowledged. 


them. The quality of the fits is comparable with the two-parameter version of the 
GCM discussed in Section [ 2 T| Fig. [T5] shows an IBM -1 fit as an example, which is 
to be compared with the experiment in Fig. Like the GCM, the standard IBM 
encompasses only the limit of 7 -instability. To stabilize triaxility, the IBM Hamil¬ 
tonian has been complemented by a third order term in Q^P^The extension is well 
exposed in Ref. ,1^ which discusses the staggering of the 7 band in ^^^di 2 ,ii 4 j^^ ^ 

signature for triaxiality. Since the appearance of triaxiality in the IBM framework 
is analog to the one in the framework of the BH discussed in Sect. | 2 . 2 [ it will not 
be further reviewed. 

5.1. Boson Number Limit 

The IBM has been conceived as an approximation to the Shell Model. The config¬ 
uration space of the valence nucleons between two closed shells is truncated to the 
subspace of pairs coupled to spin zero and two, which are then mapped to the space 
of the s and d bosons. The number Nb of such pairs is taken one half of the number 
of valence particles below the middle of the open shell and one half of the valence 
holes above the middle. The finite boson number is considered to be the major 
difference between the IBM and the phenomenological versions of the BH, which 
can also be cast into algebraic form (see the profound discussion in the textbook by 
Rowe and WoocP). The consequences of the boson number limit have not been well 
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Fig. 16. Energies (left) and B{E2,1 -^1 — 2) values (right) of The GCM calculation, 

denote d by ACM, is the same as the one shown in Fig. and the IBM calculation is the same as 
in Fig. |l5| Data from Ref.I^ 


studied because the IBM is usually applied to low-spin states, which correspond 
to bosons numbers far below the limit. The recent measurement of the life times 
of the yrast states in ^^^Pd in Ref.^^ provides such test. As illustrated by Fig. 16 


the yrast levels form a regular sequence of collective excitations with increasing 
B{E2^ 1^1 — 2) values up to / = 14. The data are compared with the GCM and 
IBM calculations that give the spectra shown in Figs. |5] and [Ts) According to the 
IBM counting rule, ^g^Pdse has a boson number of Nb = 5, four valence proton 
holes and six valence neutrons with respect to Z = A" = 50. The finite number of 
bosons limits the regular yrast sequence at / = 10 , which is the maximum that 
can be generated by five d bosons. The B{E2) values decrease toward the limit 
where they are zero. These consequences of the finite boson number are in clear 
contradiction with experiment. The GCM calculation, which does not assume a 
boson limit reproduces the experiment quite well. Figs. and demonstrate that 
both the GCM and the IBM reproduce the low-spin part of the spectrum, where 
the boson cut-off is of minor importance, with comparable accuracy. The example 
shows that the collective angular momentum cannot originate from the spherical 
valence particles and holes coupled to spin zero and two only. Testing the conse¬ 
quences of the finite boson number assumption of the IBM for other cases would 
be interesting. 


5.2. Mean field mapping 

IBM assumes that the s and d bosons represent valence nucleon pairs in spherical 
orbitals that are coupled to spin 0 or 2. A microscopic derivation of the IBM param¬ 
eters starting from this concept has not been succeeded for nuclei located far in the 
open shell. An alternative approach by Nomura and collaborator^^ and Bentley 
and FrauendorP^ has provided encouraging results. The IBM parameters are deter¬ 
mined by adjusting the IBM potential energy surface (PES) Ejbm{Pb^7b) given 
by Eq. <§ to the PES EmfiPF^lF) calculated by means of constraint mean field 
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theory. For simplicity, the version of Ref.^^ is presented, which determines the pa¬ 
rameters of the IBM -1 Hamiltonian ( [ 2 ^ . The mean field PES is calculated by means 
of the micro-macro approach. The mapping contains four free parameters, the en¬ 
ergy scale ce^ the deformation scale /3b = cb/3f^ and the IBM parameters y and 
The deformation scale is fixed by the requiring that the location of the min¬ 
imum of Eibm{/3b^Jb) agrees with the location of the minimum of Emfi/SF^Jr)- 
The energy scale ce is fixed by the energy of the 2^ state, which is calculated by 


means of the cranking model (see Section 6.1) or taken as the experimental value. 
The IBM parameters y and ( are found by fitting Ejbm{Pb^7b) to EmfiPF^jF) 
in the region E < IMeV, which is essential for the low-energy collective states, and 
assuming jb = 7f- Ref.calculated the IBM parameters and spectral properties 
for the nuclides with (Z, N) = (36, 40-50), (42, 56-66), (46, 56-70), (48, 60-70), 
(64, 88-98), ( 66 , 90-100), ( 68 , 88-100). The method describes the spectral structure 
of the nuclei amenable to the IBM parametrization fairly well. If c^; is fixed by the 
experimental energy of the 2^ state the spectrum tends to be too stretched. Using 
the cranking value results gives a scale close to experiment. 

Nomura and collaborators used the IBM-2 version. They introduced an addi¬ 
tional term oc /(/ + 1 ) the coefficient of which is determined by a cranking calcula¬ 
tion. In fitting Eibm{Pbi1b)^ they assumed that the neutron and proton parame¬ 
ters are equal. Calculations mapping Emfi^F^lF) generated from two versions of 
a constrained Skyrme EDF have been carried out for the nuclides with (Z, N) = 
(62, 84-94), (56, 54-66), (54, 54-66), (44, 54-80), (46, 54-80), (74-76, 130-140) in 
Ref.pni(62, 84-96), (92, 144-146) in Ref. ,>^ (64-66, 84-96) in Ref.^^ Calculations us- 
ing the constraint HEB applied to the Gogny interaction have been carried out for 
the nuclides with (Z, N) = (78, 102-120) in Ref. ,1^ (74-76, 114-120) in Ref.,l^(70- 
78, 110-122) in Ref.^^ The experimental energy scale and the spectral structure is 
well reproduced for the nuclei amenable to the IBM phenomenology. Like in Ref.,1^ 
the Oj states are frequently predicted too high. This concerns not only the very 
low-lying Oj states that are excluded as ’’intruder states” from the practiced IBM 
phenomenology, but also many of the Oj states in well deformed nuclei. The Oj 
states in X(5) nuclei are reproduced best, because they have the character of a soft 
anharmonic vibration in the p degree of freedom. In Ref. ,1^ Nomura et al compared 
for i92-i96p^ results obtained by mapping EmfiPF^F) from RMF to IBM-2 
with 5DBH-RMF calculations based on the same RMF. The IBM-2 generates a 
quasi 7 band with the pronounced even/-low staggering of 7 -soft nuclei whereas 
5DBH-RMF, consistent with experiment, gives very small staggering signaling the 


stabilization of the triaxial shape (cf. Sec. 2 . 2 ). One may expect that calculating 
the parameters of an IBM Hamiltonian augmented by a third order term by means 
of the mapping procedure will fix this deficiency. 

The success of the mapping technique appears puzzling from the point of view 
of the ATDF approach, because the values of the mass parameters are determined 
by the the potential energy surface. The relation, which has not been made explicit 
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yet, originates from the algebraic structure of the IBM. However, Fig. does not 
reveal an obvious correlation between and H, which is an example for the 
general situation. 


6. Non-adiabatic approaches 

The various approaches in Sections i-m presume that the quadrupole mode is 
decoupled from the quasiparticle excitations. Fig. schematically illustrates that 
this presumption does not hold far. The Oj member of the two-phonon triplet is 
usually close to the lowest 0 + two-quasiparticle excitation and some coupling is 
expected,. (The situation for the (3 vibration in well-deformed nuclei is analogous. ) 
The vicinity of the quasiparticle excitations is the reason why the collective models 
perform poorest for the 0+ excitations. The low-/ members of the higher phonon 
multiplets become progressively mixed with the quasiparticle excitations, as seen in 
the E 2 strength distribution of the Shell Model calculation on the right side of Fig. 

Very soon soon they are ” dissolved in the sea of quasiparticle excitation”, which 
means they cannot be described by a coherent collective wave function. According 
to BH phenomenology, the 0 J states of well deformed nuclei are expected to be a /3 
vibration around the axial equilibrium shape. However, the careful analysis by Gar- 
rettP^ demonstrated that the properties of the 0 J excitations deviate qualitatively 
from the characteristics of a /3 vibration. Only the the very collective Oj states 
in transition point nuclei come close to the predictions of the phenomenology or, 
alternatively, can be interpreted in terms of shape coexistence. 

For the yrast states of the multiplets the sea of quasiparticle excitations is 
approached more slowly, because the level density increases as a function of the 
distance from the yrast line. In well-deformed nuclei, the low part of the yrast line 
is just the rotational band built on the deformed ground state. Around / = 10 this 
ground (g) band encounters the s band, which is built on a non-collective excitation 
composed of two high-j quasiparticles which align their spin. This interplay between 
the collective rotational degree of freedom and the quasiparticle degrees of freedom 
can described in considerably detail in the framework of the rotating mean field 
approaches. The deformation of the nucleus is considered as static. Depending on 
the specific version it is optimized for each angular momentum (’’Total Routhian 
Surfaces” TRS, ’’Cranked Nilsson Strutinsky” approach CNS) or kept constant 
(’’Cranked Shell Model” CSM). The coupling between the collective rotation and 
the quasiparticle motion is taken into account in a non-adiabatic way. The present 
article will not address this field (see FrauendorP^ and Satula and WysP^for recent 
reviews). Two approaches for near-yrast states will be reviewed: the tidal wave 
concept that allows one to describe non-adiabatic coupling between the collective 
quadrupole mode and the quasiparticle modes in the vibrational and transitional 
regime and the Triaxial Projected Shell Model that allows one to take into account 
the collective 7 mode. 
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Fig. 17. Left panel: The potential 7e) along the path of the deformation minimum. 

Right panel: Experimental moments of inertia of the ground bands, of the Nd isotopes, where 
ca = {E{I) - E{I - 2))/2 and J = I. From Ref.I^ni 


6.1. Tidal Waves 

Frauendorf, Gu and Sun introduced the Tidal Wave concept in Ref!^ and Ref.,1^ 
which contains complimentary material. 


6.1.1. Phenomenology 


Consider the classical BH ( 12 ) with phenomenological mass parameters Bgg = 
= Bi = y/b/2D. Uniform rotation about the medium (m-) axis, which has the 
maximal moment of inertia, has the lowest energy for a given angular momentum, 
i. e. it corresponds to the yrast state when quantized. The choice of the interval 
27r/3 < 7 < 7 r /3 makes the m-axis the 3-axis of quantization. The location of the 
surface in spherical coordinates is given by 


R{t) = Ro[l V 2/3 sin7cos(20 — 2ujt)Y22{d, 0 = 0). 


(28) 


The deformation parameters (3 and 7 do not depend on time, because any time 
time dependence involves additional kinetic energy. As discussed in detail in Ref.,1^ 
their values are given by minimizing the energy 


^(A7) = 


J" 


2J(A7) 

Let us discuss some cases in more detail. 

U 
2 


+ !/(/?, 7 ), J = 45/32 sin 7 . 


-Harmonic vibrator:(V in Fig. 17) V = ^(3^ 
Minimizing the energy one finds 


7e 


- 0^ 
2 ’ 


J 


2v^ 


, J = AB0l = 


2J 




(29) 


(30) 


( 31 ) 
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J 


Fig. 18. Experimental moments of inertia of the ground bands, of the Cd, Pd, and Ru isotopes, 
where cj = {E{I) - E{I - 2))/2 and J = 1. 


The wave travels with an angular velocity uj being one half of the oscillator fre¬ 
quency Q. The angular momentum is generated by only increasing the deformation 
which increases the moment of inertia, which is a linear function of J. This 
mode has been called ’’tidal wave”, because it has wave character: the energy and 
angular momentum increase with the wave amplitude while the frequency stays 
constant. 

-Axial rotor (R in Fig. [Tt])): The potential has a minimum at ydo, where 
V ^ C(/3 — /do)^/2 and is stiff around je = tt/S. 

The minimization gives 

n t2 o t2 

/3e = /3o(l + + 0(j% J = m + + 0{T)). (32) 

The moment of inertia is a slowly increasing quadratic function of J. Energy and 
angular momentum increase (mainly) due to growth of the angular velocity uj. 
-Transitional nucleus (T in Fig. [T^): The yrast energies are well approximated 
(J = /) by the expression 

J 

= E{J) = co{J)dJ. 


fj — ©0 + ©iT, cj 


(33) 
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Assuming that the energies result from minimization of the energy (29), the poten¬ 
tial V{Peile) is given in parametric form by 


PeiJ) sin7e(J) 


00 + ©1J 


4B 


V{J)= [ u{J)dJ-u{J)J/2, 

Jo 


(34) 


which has an intermediate form T in Fig. 17 As expected for a transitional nucleus, 
the angular momentum is gained by increasing both J and uj. From a vibrational 
perspective, the increase of uj reflects the anharmonicity of the motion, from the 
rotational perspective, the increase of J' reflects the softness of the rotor. 

The chain of Nd - isotopes in Fig. [^displays the transition from a vibrational 
(V - small ©o) to a rotational (R - large ©o) yrast sequence. The value of 7e cannot 
be inferred from the yrast sequence only. In the case of the Nd isotopes the high 
energy of the 7 band points to a near-axial value close to 27r/3 (X(5) type). For 
A" = 84, the vibrational sequence is short and somewhat irregular, which indicates 
that the underpinning fermonic structure seems through the weakly collective mode. 
For A = 86, the increased collectivity smoothes out most of the irregularity, which 
disappears for A = 88. The Cd chain in Fig. has small values of ©0 as expected 
for nuclei close to the vibrational limit. The quasiparticle excitations are early 
encountered. The isotonic Ru chain has transitional character with 5 MeV< ©1 <15 
MeV. The low ^(2 J) energy and the pronounced even-/-up staggering of the 7 band 
point to 7e ~ 7r/2 (E(5) type). It is obvious that the encounter of the quasiparticle 
excitations does not correlate with the pair counting rule of IBM. 

Semiclassically, B{E2,I I — 2) oc Qf, because the transition quadrupole 
moment Qt = Q 2 ^ Pe sin 7e according to Eq. 0. Therefore the ratio B{E2)/J 
does not depend on J(= I). This correlation is well known as Grodzin’s rule for the 
2^ states. It also applies with good accuracy for the I > 2 yrast states in nuclei with 
44 < Z < 48 as far as the lifetimes of the states have been measured. Eig.[T^ shows 
^^^Pd as an example. Other regions have not been systematically investigated in 
this respect, though it holds for Gd and ^^^Pt. 

Ayangeakaa et alW and Macchiavelli et alW^ interpreted the classical tidal wave 
as a condensate of d bosons. Accordingly, up to seven bosons are observed, which 
align their angular momenta. As discussed in Section |5.H the structure of these 
d bosons must be more complicated than pairs of valence particles in a spherical 
potential coupled to angular momentum of two. If the bosons were free, the function 
J{I) would be a straight line out of the coordinate origin (EB in Eig. 19 and V 
in Eig. prl). The displacement by ©0 was attributed to an interaction between 


the bosons (IB EB in Eig. 19) that is quadratic in the boson number. Erom the 
systematics in the region 44 < Z < 48 a correlation was identified: the smaller the 
boson energy the larger is the boson interaction. 
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Spin I (ft) 



Fig. 19. Left part: Experimental moments of inertia of the ground band of ^®^Pd, where 
w = {E{I) — E{I — 2))/2 and J = I. Right part: Experimental B(E2^ I -^1 — 2) values of ^®^Pd. 
The blue lines WS-SCTAC show the calculations by means of the cranking model.^^ From Ref.ini 


6.1.2. Semiclassical microscopic calculations. 

Semiclassically, the yrast states of the quadrupole mode are represented by a tidal 
wave that travels with constant angular velocity uj over the nuclear surface. This 
has allowed Frauendorf, Gu,and Suii221l£il to calculate its properties by means of the 
SC-TAC model (for details se^^. The potential 'F(/3,7) in Eq. (29) is determined 


by the Micro-Macro method described in Section [3^ It is based on the mean field 


Hamiltonian h^BCS^ which has the structure of hpQQ given by Eq. (16), combining 
the Nilsson deformed potential with the monopole pair field. The rotational energy 


in Eq. (29) is obtained by the cranking procedure 




= E'{uj,/3,j)+ojJ - E'{0,/3,j), {uj,/3,j\j^\oj,/3,'y) = J (35) 

[hNBCsiPn) -(^jx] | w ,/ 3 , 7 ) = £''( w ,/ 3 , 7 )| w ,/?, 7 ). (36) 


Eor each deformation grid point /3,7, the frequency uj{J) is adjusted to obey the an¬ 


gular momentum constraint (35). The mean field state jcj, 7) is found by the solv¬ 


ing the the quasiparticle eigenvalue problem (36) exactly. This is a crucial difference 
to the ATDME, which takes the cranking term 00 jx into account by perturbation 
theory, resulting in the IB expression for the moment of inertia. The equilibrium 
values /de, 7e are found by minimizing E( J, 7). The transition quadrupole moment 
is obtained as the expectation value (J, /de, 7e|Q2|'/,7e) 

Erauendorf, Gu,and Suil221l£ll calculated the energies of the yrast states and 
the B{E2) of the intra band transitions up to spin / = 16 for the nuclides with 
Z = 44—48, N = 65—66. The ^-factors for the same states were calculated in Ref.,1^ 


where the experimental deviations from Z/A could be reproduced. Eigs. 19 and[^ 
exemplify the accuracy of the parameter-free calculations. In particular the change 
of the yrast states from the purely collective tidal wave (g band) to the configuration 
with two rotational aligned hi 1/2quasiparticles (s band) is reproduced in detail. In 
























June 23, 2015 0:26 WSPC/INSTRUCTION FILE 


ws-ijmpe-Frauendorf 


32 S. Frauendorf 



o[MeV] 


"a 



I 


Fig. 20. Left part: Experimental moments of inertia of the ground band and Oj band (black dots) 
of ^^°Cd, where uj = (E(/) — E{I — 2))/2 and J = I compared with the calculations by means 
of the Cranking modePS (black circles). Right part: Experimental B(E2, / -^1 — 2) values of 
^^°Cd compared with the calculations (black circles without red circle with quantal correction, 
see Ref.I^). From Ref.^^ 


the case of (Fig. 19) the collective g band can be followed up to / = 14, where 

it is at higher energy than the s band. This is a consequence of almost no mixing 
of the two configurations. In the case of ^^^Cd (Fig. 20) the collective g band is 
crossed by the s band earlier and the two bands interact stronger. The two aligned 
hii /2 quasiparticles in the s band reduce the deformation but stabilizes it such that 
the sequence becomes more rotational. The rotational sequence based on the low- 
lying Oj state is also well reproduced. In the calculations it is a two-quasi proton 
excitation with larger deformation. Its calculated energy is substantially larger than 
in experiment, which points to missing quadrupole correlations of vibrational type. 
More systematic calculations for other regions would be of interest. The method 
applies to odd-A and odd-odd nuclei without any further sophistication. 


6.2. Triaxial Projected Shell Model 

Sheikh and HargP^ introduced the Triaxial Projected Shell Model (TPSM) to de¬ 
scribe the quadrupole mode in triaxial nuclei. The TPSM is based on the earlier 
axial version of the PSM, which has widely used for interpreting high spin exper¬ 
iments. Since this work has been reviewed in Ref.,1^ only the more recent work 
on triaxial nuclei will be covered here. The TPSM starts with a superposition of 
a set of quasiparticle configurations |/i:) in a triaxial Nilsson potential, which are 
projected by means of Plfx good total angular momentum I and its projection 
K on one of the principal axes of the potential. 
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Fig. 21. The spectr a of Nd and ^^^Ce compared with TPSM calculations. Reproduced with 
permission from Ref.^^ 


\IMa) — 

kK 

k'K' k'K' 

'HkK,k,'K' = {i^\HpqqPxk'\'^') 5 -^k,k,k,'K' = {e\Pkk'\'^') • 


(37) 

(38) 

(39) 


The amplitudes are found by diagonalizing the Pairing plus Quadrupole Hamil¬ 


tonian HpQQ Eq. (15) in the non-orthogonal basis. They represent the collective 
quadruple mode of the triaxial rotor, which can be excited into multi-quasiparticle 
configurations. The axial deformation parameter p is taken either from the mea¬ 
sured H(E2, 2+^0+) values or calculated in the framework of some mean field 
approach. It determines the coupling constant p of the QQ-interaction. The triaxi- 
ality parameter 7 is adjusted to reproduce the energy of the 2^ state (head of the 
quasi-7 band). Values close to 30*^ are found for all cases. 

The coupling of the quadrupole mode to the quasiparticle excitations is a new 
quality. The strong coupling between the collective and quasiparticle degrees of 
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freedom leads to new physics, as e.g. the appearance of chirality in ^^^Cs, which 
can be described in the framework of the TPSM!^ The capability of TPSM to 
incorporate the coupling to the quasiparticle degrees of freedom makes it an efficient 
tool to describe the structure of odd-A and odd-odd nuclei on a microscopic basis. 
Examples are the studies of Ref.l^and ^^^Nb in Ref.^^ 

Sheikh et al studied in detail the properties of the 7 band in the isotope chain 
156-172^^ in Refs.flS^K^The calculations very well reproduce both energies and in¬ 
traband B{E2) values. The work of Boutachkov et demonstrates that TPSM 
also reproduces the transition probabilities between the 7 and ground state band 
rather well. The calculations include two-quasineutron and two-quasiproton excita¬ 
tions and the combinations thereof. This makes it possible to describe the alignment 
of high-j quasiparticles on top of a 7 vibration, which is out of reach of other ap¬ 
proaches. In the TPSM framework, the low-lying K = 3 bands observed in these 
nuclei are interpreted as built on triaxially deformed two-quasiparticle states.^^^ 
TPSM is predestinate for describing triaxial nuclei. RefP^ investigated 
i80-i90p^ and suggested that the low-lying Oj state is a two-quasiproton excitation. 
This is in contrast to earlier interpretations which assign an axial shape to it that 
differs from the ground state shape. The rigid-rotor like spectrum of ^^^Ru is very 
well accounted for by the calculations in Refill The band structures of 
and i32-i38]N^q accurately reproduced in the TPSM calculations of Refill Fig. 
2l\ shows how the calculations describe both the collective ground and 7 bands 
and the high spin bands built on configurations with a rotationally aligned proton 
or neutron pair. The study established a collective excitation of the 7-vibration 
type built on such two-quasineutron configuration, which explains the negative g 
-factor of the 10+ state. Ref-EnH addressed the even-odd staggering of the 7 - band 
in and ^^“^^Se. As discussed in Sect. 2.2, the phase of the staggering pa¬ 


rameter S{I) indicates if the nucleus is 7-soft or rigid in the framework of the BH 
phenomenology. The TPSM calculation reproduces the surprising observation that 
S{I) has the low-odd-/ character of 7-stiff nuclei for ^^Ge while all neighbors have 
the low-odd-/ character of 7-soft nuclei. 

The reason for the success of the TPSM in accounting for the details of the 
quasi-7 bands remains to be understood. The quasi-7 vibration can be understood 
as a wave traveling the axial-symmetric surface of the nucleus, very much like the 


tidal wave the surface of a spherical nucleus (cf. Sec. 6.1). This corresponds to a 
static triaxial deformation in the co-rotating frame of reference. However, it seems 
counterintuitive that all considered cases, encompassing the relatively high-lying 
7 bands in the well deformed Er isotopes and the low-lying quasi-7 bands in 7- 
soft transitional nuclei with both phases of the staggering parameter 5'(/), are well 
described by a deformed mean field with 7 ^ 30^. The properties of the quasi-7 
bands do not correlate with the mean field value of 7 in any obvious way, which 
is in stark contrast to the BH phenomenology. Moreover, a certain admixtures of 
two-quasiparticle configurations is needed to describe the experiment. 
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The authors show the energies of the basis states before mix¬ 

ing, which they call ’’band diagrams”. However the insight provided by these graphs 
is limited, because the properties of interest often emerge as the consequence of the 
strong mixing a number of such basis states. A more sophisticated analysis of the 
TPSM wave functions seems desirable to understand the physics behind the TPSM 
results. 

7. Assessment and challenges 

Basing the phenomenology on the BH with irrotational flow type mass parameters 
and a forth oder potential (GCM) provides a classification scheme that qualitatively 
accounts for many features of the quadrupole mode in terms of three parameters. 
The model provides an illuminating interpretation in terms of the nuclear shape 
for data in the whole range between the text book limits of harmonic vibration, 
rigid rotation (axial and triaxial), and 7 instability. The same holds for the IBM 
phenomenology, which has become popular because of the availability and simple 
use of the IB A codes. A similar handy GCM code for public use would be very 
welcome. Both models are used to describe the properties of the experimental 2^, 
4^, 6^, 2^, 3^, 4^, Oj, 4^ states, which they reproduce with comparable quality. 
The often claimed foundation of the IBM on the microscopic Shell Model cannot 
be substantiated. The boson counting rule, which follows from this assertion, has 
not been systematically tested, and it fails for ^^^Pd. The derivation of the IBM 
parameters using truncation schemes of the Shell Model configuration space have 
not succeeded so far. Instead mapping of the potential energy surface calculated 
from the deformed mean field to the surface generated by the coherent state repre¬ 
sentation of IBM results in parameters that well account for the experiment. 

The microscopic form of the BH is obtained by applying the adiabatic approxi¬ 
mation to the TDMF theory. There are several versions that start from the various 
mean field approaches used in practice at present, which are carefully tuned to 
experimental quantities other than the collective quadrupole excitations. The po¬ 
tential energy of the BH can be considered as reliable as the mean field approaches 
reproduce the properties of the single particle levels near the Fermi surface. In this 
respect, the modern versions based on the Skyrme Energy Density Functional, the 
Relativistic Mean Field, and the Gogny Interaction are not superior to the older 
ones that start from a Modified Oscillator or Woods Saxon potential, which are 
directly adjusted to experimental properties of the single particle levels. Not sur¬ 
prisingly, the accuracy of the different approaches is comparable, because it cannot 
be better than the underpinning mean field theory. All versions but one use the 
IB form for the mass parameters, which gives too small values for the mean field 
treatment of the pair correlations, resulting in a too dilute excitation spectrum. 
The discrepancy is removed by either simply scaling the spectrum or reduction of 
the pair correlations. The models perform well for the yrast and yrare states with 
I < 6 with considerable predictive power. The performance is poorest for the ex- 
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cited 0+ states, which are expected to couple strongest to the quasiparticle states. 
Nevertheless, low-lying 0+ states with shape coexistence, vibrational character or 
intermediate character are quite well accounted for. The pronounced deformation 
dependence of the mass parameters is essential for establishing these structures. 
Calculating by means of the microscopic BH the spectroscopic properties of the 
quadrupole excitations across the mass table is feasible with modern computers. 
Such a project has been only carried out for the version based on the Gogny inter¬ 
action, the results of which is accessible online. To create a similar resource for the 
other versions of the microscopic BH would be very useful, both for interpreting 
experiments and for a more systematic assessment of the model performance. 

Taking into account the time-odd parts of the quadruple pair field by means of 
the LQRPA method increases the mass parameters from the IB to the TV values, 
which give the correct energy scale. The LQRPA is computationally demanding, and 
it has only been applied for the simple PQQ interaction for this reason. Developing 
tractable methods to calculate the TV mass parameters for the other types of 
interactions or EDF’s is a challenge. 

The Generator Coordinate Method (without GOA) applied to all five 
quadrupole degrees of freedom is computationally very extensive. For this reason 
it has been restricted to very light nuclei. The method is particularly well suited 
for parallel computing, which may lead to progress with the further development of 
super computers. For the heaviest nucleus studied, ^^Kr, it gives similar result as 
the microscopic BH derived from the same RMF EDF. So far the collective states 
have been generated from the ground state configurations on the quadrupole coor¬ 
dinate grid, which implicitly invokes the adiabatic approximation. The coupling to 
quasiparticle excitation can be incorporated in a straight forward way, such going 
beyond the adiabatic regime. To make this tractable is certainly a challenge. 

A wide stride into this direction is the Triaxial Projected Shell Model, which 
assumes a fixed deformation, projects on good angular momentum, and takes cou¬ 
pling to up to four-quasiparticle configurations into account. Using the simple PQQ 
interaction makes the calculations practical. The method includes an element of 
phenomenology, because the shape of the mean field is partially adjusted to re¬ 
produce the experiment. The coupling between the collective quadrupole and the 
quasiparticle excitations is remarkable well reproduced for the states in the yrast 
region. This is also true for odd-A and odd-odd nuclei, to which the method ap¬ 
plies without any complications. Often the results are obtained as a heavily mixed 
number of components in a non-orthogonal basis. It is a challenge to find inter¬ 
pretations of such results in simple terms in oder to elucidate the physics behind 
them. Treating the deformation parameters as generator coordinates appears an 
attractive future development. The simplicity of the PQQ interaction will keep the 
computational requirements as low as possible, which facilitates the exploration 
of the new territory. The not-yet-understood observation that the model describes 
the various modes of the quasi-7 band (7-soft vs. 7-rigid) assuming one static mean 



June 23, 2015 0:26 WSPC/INSTRUCTION FILE ws-ijmpe-Frauendorf 


The Low-Energy Quadrupole Mode of Nuclei 37 


field triaxiality of 7 ^ 30^ may considerably simplify the calculations. 

The Tidal Wave concept provides a simple semiclassical description for the yrast 
states. The interplay between the collective quadrupole mode and the quasiparticles 
is described by quasiparticle configurations in a uniformly rotating potential. The 
method works for the whole range of vibrational to rotational nuclei, for which it 
is well tested and used as the standard classification of the rotational bands. The 
method very well reproduces the energies of the yrast levels and the B{E2) values of 
the connecting transition for the Ru, Pd, and Cd isotope chains. It can be directly 
applied to odd-A and odd-odd nuclei. The semiclassical approximation is not quite 
accurate for the low- spin states / < 6. Projecting on good angular momentum 
should remove this deficiency. 

Exploring the region of strong coupling between the collective and quasiparticle 
degrees of freedom will require a paradigm change. Some large-dimension diagonal- 
ization has to be used to account for the rapid increase of the level density. The 
nucleus is a partially chaotic system. Only close to the yrast line one can expect 
that the theory reproduces the experiment state by state, and the traditional com¬ 
parison of individual energies and transition probabilities is appropriate. Higher up 
one has to take into account a certain degree of randomness and study the proper¬ 
ties of groups of levels. It is a challenge to formulate new concepts that identify the 
correlations due to the fragmented quadrupole mode as long as it survives. This 
concerns also standard Shell Model calculations like the example in Fig. It is not 
clear from simple plot of the B{E2) values of the connecting transitions whether 
the low-spin members of the three- and higher phonon multiplets survive in a frag¬ 
mented form or if the quadrupole collectivity is quenched. The availability of the 
full information about the wave function allows theory to address this question in 
novel ways. 
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